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PROCEDE DE DEMODULATION AVEUGLE AUX ORDRES SUPERIEURS 
D'UN EMETTEUR DE FORME D'ONDE LINEAIRE 

5 Uobjet de Tinvention concerne un precede de demodulation 

aveugle de signaiDc 6m\s par plusieurs ^metteurs et reQus par un r^eau 

compost d'au molns un capteur. 

II s'applique par exempie pour un r^eau d'antennes dans un 

contexte Slectromagn^tique. 
10 Uobjet de ('invention concerne notamment la demodulation de 

signaux, c'est-§-dire Textraction des symboles {aj emis par un dmetteur 

moduli lindairenrrent. 

La figure 1 reprSsente un systSme de traitement d'antennes 

15 comportant plusieurs 6metteurs Ei et d'un systeme T de traitement 
d'antennes comportant plusieurs antennes Rl recevant des sources radio- 
electriques avec des angles d'incidence diffgrents. Les angles d'incidences 
des sources ou emetteurs peuvent etre paramStr^s soit en 1D avec Tazimut 
Om soit en 2D avec Tangle d'azimut Om et Tangle d'elevation Am. 

20 La figure 3 schematise un principe de modulation et demodulation 

des symboles { ak } Smis par un dmetteur. Le signal se propage au travers 
d'un canal d multi-trajets. Uemetteur emet le symbole ak d Tinslant k.T, oCi T 
est la pdriode symbole. La demodulation consiste d estimer et d detector les 
symboles pour obtenir en sortie du demodulateur les symboles estimSs Sk- 

25 Sur cette figure, le train de symboles {ak} est filtre lineairement d I'emission 
par un fiKre d'emission H appele egalement filtre de mise en forme ho(t). 

Dans la suite de la description, on definit sous Texpression 
« demodulation aveugle », des techniques qui n'utilisent pas d'information a 
priori sur le signal emis : filtre de mise en forme, sequence d'apprentissage, 

30 etc.. 
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Les dix demidres annees ont vu le developpement des techniques 
de demodulation aveugle SIMO, abr§g6 de entree unique sorties multiples, 
(en abrige anglo-saxon Single Input Multiple Output) dites d sous-espace 
utilisant les statlstlques d'ordre 2, telles que d^crites dans la r^rence [7]. 
5 Ces algorithmes presentent toutefois rinconv^nient de ne pas Stre robustes 
ni d une sous-estimation ni d une sur-estlmatlon de I'ordre du canal de 
propagation : statement temporel dependant des multi-trajets et du filtre de 
mise en forme. Pour contourner ce probldme, II a 6t6 propose une technique 
de prediction linSaine ddoite dans la rdf^rence [11] qui pr^sente comma 
10 inconvenient d'etre moins performanta lorsque la longeur du canal est 
connue. Pour ameiiorer les techniques d sous-espaoe, la methode decrite 
dans [16] propose une technique parametrique qui necessite 
malheureusement la connaissance du filtre de mise en fonne. 

Dans la reference [13], les auteurs proposent une technique § 

15 base de covariance matching, ayant notamment I'inconvenient d'etre tr6s 
difficile a mettre en oeuvre. C'est ainsi qu'il a et6 developpe une technique 
sous-optimaie decrite dans la leterence [12] plus facile d mettre en csuvre en 
minimisant un critere de vraisemblance et supposant le caractere gaussien 
des symboles. Cette hypothdse n'est pas verifies pour les modulations 

20 lineairas couramment utilisees telles que les PSK (Phase Shift Keying) ou les 
QAM (abieviation anglo-saxonne de Quadrature Amplitude Modulation). 

II est aussi connu dans les methodes CMA (en abi^e anglo- 
saxon Constant Modulus Algoritiim) d'utiliser une approche spatio-temporelle 
decrite par exemple dans la reference {6}. Ceite famille de methodes 

25 presente toutefois i'inconvenient de n'etre adaptee qu'd une dasse 
parSculiere de modulations telles que les PSK qui sort d module constant 
Cette methode est rterafive et a done I'inconvenient de devoir §tre 
correctement initialisee. Pour finir, les meUiodes CMA presentent le 
desavanfage de converger moins vfte que la methode d sous-espaoe 

30 enoncee precedemment D'autre part, la reference [20] decrit une methode a 
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sous-espace exploitant les statistiques d'ordres superieurs pour 
ridentification de canaux a RSponse impulsionnelle Finie (FIR) et d phase 
non minimum. 

5 L'objet de la presente Invention oonoeme un proo6dd bas6 

notamment sur des techniques de separation de sources en aveugle 
connues de i'Homme du metier et ddcrites par exempie dans les references 
[4] [5] [15] [19] supposant que les symboles emis sent statistiquement 
Independents. Pour cela ie precede construit une observation spatio- 
10 temporede dont les sources meiangees sent des trains de symboles de 
remetteur. Cheque train de symboles est par exempie Ie mSme train de 
symboles decaie d'un nombre entier de periode symbole T. 

L'invention oonceme un precede de demodulation aveugle d'une 
source ou emetteur de forme d'onde lineaire dans un systdme comportant 
15 une ou plusfeurs sources et un rSseau de capteurs et un canal de 
propagation caracterise en ce qu'il comporte au moins les etapes suivantes : 

• determiner Ie temps symbole T et on echantillonne a Te tel que T=ITe 
(I entier), 

• a partir des observations x(kTe), construire une observation spalio- 
20 temporelle z(t) dont les sources meiangees sent des trains de 

symbole de Pemetteur , 

• appllquer une methode de type ICA sur Ie vecteur d'observation z(t) 
pour es&mer les U trains de symboles { am-i } assodes aux vecteurs 
decanal 

25 • ordonner les U sorBes (am. j, j ) dans ie meme ordre que les entrees 
(am^ , afin d'obtenir les vecteurs de canal de propagation 

• determiner la phase aimax associee aux sorties. 



30 Le precede selon I'invention ofFre notamment les avantages 

suivants : 
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• II ne lait aucune hypothese sur les constellations de symboles 
contrairement aux methodes decrites dans Tart anterieur, 

• II ne ndcesstte pas la connaissance du filtre de mise en forme. 

• Le module des symboles n'est pas suppose constant, 

5 • II est robuste d une surestimation de la longueur du canal, 

• 11 permet de traitor le cas des canaux de propagation avec des trajets 
corr§leS| 

• II est direct et simple d mettre en oeuvre sans 6tape de recoupement 
des trajets corrSles. 

10 

D'autres caractSrlstlques et avantages de Tobjet de la prdsente 
Invention apparaTtront mieux a la lecture de la description qui suit donn^e a 
titre illustratif et nullement limitatif a la lecture des figures annexdes qui 
repr^sentent : 
15 • La figure 1 un exemple d'architecture, 

• La figure 2 les angles d'incidence des sources, 

• La figure 3 le processus de la modulation lin^aire et demodulation 
d'un train de symboles, 

• La figure 4 le schema d'un emetteur a modulation lineaire, 

20 •La figure 5 un resume du principe general mis en oeuvre dans 
rinvention, 

• La figure 6 la representation d'une constellation, 

• La figure 7 un premier exemple de mise en ceuvre du procdde oD le 
signal est regu en bande de base, 

25 •La figure 8 un deuxi^e exemple oD le signal est regu en bande de 
base et les multi-trajets sent d6corr6l6s, 

• La figure 9 un trolsldme exemple oD le signal est regu en bande de 
base et les multi-trajets sent ddoorr6l6s par groupe. 
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Afin de mieux laire comprendre le proced^ selon i'invention, la 
description qui suit conceme un proc^d^ de demodulation aveugle aux 
ordres superieurs d'un dmetteur de forme d'onde lin^aire dans un reseau 
ayant une stmcture telle que celle d^rite d la figure 1, par exemple. 

5 Avant d'expliciter les dtepes mises en osuvre par le proc^6, on 

ddcrit le module du signal utiilsd 
Moddle du signal 6mis par une source ou 6metteur 
Modulation lin^aire 

Les figures 3 et 4 montrent le processus de la modulation iin^aire 

10 d'un train de symboles {aij d la cadence T par un filtre de mise en forme 
ho(t). 

Le peigne de symboles c(t) est tout d'abord filtrd par le filtre de 
mise en forme ho(t) et ensuite transpose d la frSquenoe porteuse fg. Le filtre 
NRZ, qui est une fenStre temporelle de longueur T , tr^s souvent difini par 

15 ho(t}=ni<t-T/2), est un exemple parficulier non limitatif de filtre d'^mission. 
Dans les radiocommunications, 11 est 6ga!ement possible d'utiliser le filtre de 
Nyquist dont la transfbrmSe de Fourier ho(f)«n8(f-B/2) se rapproche d'une 
fenetre de bande B. lorsque le roll-off est nul alors ho(f)=nB(f-B/2) (le roll-off 
definit la pente du filtre en dehors de la bande B). 

20 Le signal module So(t), 6mis par I'Smetteur, s'6crit d I'instant 1k=kTe 

(Te : periode d'^chantillonnage) en fonction du peigne de symboles c(t) : 

(1) 

So(l<Te)=2 ho(iTe) C((k-i)Te) 

i 

Prenons un temps symbole T egal d un nombre entier de fois (a 
p6riode tf^chantillonnage. T=ITe et posons k=ml+j avec 0 < j < I. Pulsque 
c(t)=Sr Br 5(t-rlTe), autfement dit, comma c(t)=au pour t=ulTe et c(t) = 0 pour 
25 t^^ulTe, les seules valeurs de i pour lesquelles c((k-i)To) est non nul verifient k- 
i=ul, c'est-a-dire telles que i=ml+j-ul=nl+j ou n=m-u. Finalement, I'expresslon 
(1)devient: 
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So(mlTo+ jTe y= 2 ho(nlTe+jTe) am^, pour 0^i<l 

Le parametre U est la demi-longueur du filtre d'6mission qui s'6tale sur une 
duree de (2Lo+1)ITe. Dans le cas particuller d'un filtre d'emission NRZ, on 
obtient Lo=0. Quant au signal emis s(t), 11 v^rifie s(t)=so(t) exp027tfot) car il est 
6gal au signal so(t) transpose d la frequence fo. Dans ces conditions 
5 I'exprssslon de 8(mlTe+jTe) est d'apr^s (2) : 

s(mlTe+jTe )= i; ho(nlTe+jTe) expa2jcfo(nl+j)Ts) a^n. expG2jtfo(m-n)ITe) 

= ^ hFo(nlTe+jTe)bmH, telque 0<j<l 
ou hFo(iTe)=h(iTe) exp02icfoiTe) et bi=aj expG2jifollTe) 



Reception des signaux sur les capteurs 

Le signal s(t) (FIG.3) Smis passe d travers un canal de propagation 
avant d'§tre re^u sur un rdseau composd de N antennas. Le canal de 
10 propagation peut se mod^liser par P muKi-trajets d'inddence 8p, de retard Tp 
et d'amplitude pp (1 < p ^ P). En sortie des antennes on a le vecteur x(t) qui 
correspond d la somme d'un melange lin^aire de P mutti-trajets et d'un bruit 
suppose blanc et gausslen. Ce vecteur de dimension A6c1 a I'expression 
suivante : 

x(t)=E Ppa(Gp)s(t-Tp)+b(t)=As(t) + b(t) 

15 OU pp est I'amplitude du p*"^ trajet, b(t) est le vecteur bruit supposd 
gaussien. a(e) est la r^ponse du rdseau de capteurs d une source 
d'inddence G, A=[ a(Gi)... a(Gp)] et s(t)=[s(t-Ti) --s(t-cp)r. En notant que Xp = 
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fpT+Atp (oD (0<ATp<T=ITe) et tp est un entier) et en utilisant I'expression (3) 
dans i'6quation (4). on obtient pour le vecteur sur les antennes : 

p A, (5) 
x(mlTe+jTe) =X Pp 3(®p) •^rat" ITe+jTe-AXp) b^^,^ + b(mlTo+jTe) 

En effectuant le changement de variable suivant Up = n + tp, le vecteur roQu 
par les antennes s'exprime : 

P (6) 

x(mlTe+jTe) =£ X pp a(®p) hF0((Up- rp)ITe+jTe-ATp) b^„^ + b(mlTe+jTe) 

5 En notant rmm = min{rp} et rmax = maxfrp}, I'equation (6) peut etre reicrite de la 
manidre suivante : 

P raa±U (7) 
x(mlTe+jTe) =X ZT pp 3(®P) "^Pott"" rp)ITe+jTe-ATp) lnd[rp4J). rT>+LOl(u) 

+ b(mlTe+jTe) 

OQ lnd[r.q](u) est la fonctlon Indicatrice usuelle (lnd[r.(d(u)=1 pour i<u<p et . 
lnd[r.<d(u)=0 autrement) ddfinie sur I'ensemble des entiers relatils d valeur 
dans I'ensemble binaire {0, 1}, caract^ris^e par lnd[r^(u) = 1 si u appartient a 
10 rintervalle [r,q] et lnd[r^(u) = 0 sinon. De ce fait, en notant v(t) le vecteur de 
canal : 

p (8) 
v(ulTe+jTe) =S pp 3(®P) fV=o((U- rp)ITe+jTe-ATp) Ind(,p4j). ,p*mi(u) 

oD lr=ulTetjTe et I'expression (5) devient : 

. . (9) 



x(mlTe+jTe) = v(ulTe+jTe) ^ + b(mlTe+jTe) 



Interference entre symboles 

Le vecteur observation x(t) issu du r6seau d'antennes § I'instant 
15 t = mlTe+ITe fait, d'aprte I'^quation (9), intervenir le symbole bm mais 
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dgalement les symboles bmni ou u est un entier relatif appartenant a 
rintervalle [rmm-U. rmax+U]. ph6nom6ne qui est plus connu sous le nom 
dlnterf^rence Entre Symboles (lES). Notons U le nombre de symboles 
participant d TIES et bomons rintervalle de valeurs prises par ce dernier. 
5 D'aprds rdquation (9). si Tlntersection des intervalles [rp-U. rp+U] est non 
vide, alors on a U = I rmax - rmm I + 2Lo + 1- De ce felt, iorsque rmax = rmin. 
c'est-^d-dire Iorsque, tous les multi-trajets sent corr^lds, la borne minorante de 
U est atteinte et vaut U - 2Lo + 1- Ce cas se tradult Sgalement 
mathdmatiquement par| max{tp}-min{Tp} | < T, Par centre, si rintersection 

10 des dits intervalles est vide, et que le cas echeant, tous les Intervalles [rp-U, 
rp+Lo] sent disjoints, alors on a Lc=Px(2Lo+1)» ce qui constitue une bome 
majorante d Tensemble de valeurs susceptibles d'etre prises par U- Ce 
dernier cas de figure correspond concretement au cas de multi-trajets tous 
dScon-elSs deux k deux, ce qui matli6matiquement peut 6galement s'ecrire 

16 VI 96 j, I n - rj |>2Lo , condition obtenue des que I xi - xj | > (2Lo + 1)T, Pour 
r^sumer, la quantite U v^rifie de manlere generate Tencadrement suivant : 



Uexpression traduisant le vecteur regu par les capteurs peut alors se rddcrire 
de la maniere suivante, oD cette fois n'apparaissent que les U symboles bm-u 
d'interet : 



2Lo+1 <Lc<Px(2Lo+1) 



(10) 




(11) 



20 OD V 1 <l<Lc, et rmin-Lo< n(l) < rmm+U et oQ: 



h(t)=£ Ppa(Gp)hFo(t-Tp) 



p 



(12) 



Techniques ICA 

Le procddd lait appel d des techniques ICA bashes sur le module 
suivant donnd k titte illustfatif et nullement limitatif : 
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L (13) 

oQ Uk est un vecteur de dimension Mxl regu a I' instant k, Sir est la 
composante du signal Sk a I'instant k, nk est le vecteur de bruit et G =[ gi ... 
g J. Les mithodes ICA ont pour objectifs d'extraire les l=L composantes Sik et 
d'identifier leurs signatures g\ (La reponse vectorlelle de la source I au travers 

5 de I'observation Uk) ^ partir des observations ur. Le nombre l=L de 
composantes doit etre inferieur ou egal d la dimension M du vecteur 
d'observatlon. Les methodes des references [4] [5] et [15] utilisent les 
statistiques d'ordre 2 et 4 des observations Uk- La premiere 6tape utilise les 
statistiques d'ordre 2 des observations Uk (ces observations peuvent etre 

10 fonctions des signaux re^us sur les capteurs) pour obtenlr una nouvelle 
observation Zk telle que : 

L ^ (14) 

Zk = Wi Uk =2^ Qi sik + fik = G Sk + nk 

ou les signatures gi (1<i<L) sont orthogonales, 6 =Igi --- gj et Sk=[ Sik — 
SlkI^- La deuxieme etape consiste a identifier la base orthogonale des 6 d 
partir des statistiques d'ordre 4 des observations blanchies Zk- Dans ces 
16 conditions on peut extraire les signaux Sk en effectuant : 

Sk = 6*Zk=e*WiUk (15) 

Ou §k est Testimde des signaux Sk et oQ * est I'opdrateur de pseudo-Inversion 
d6flnipar<S*=(e"e)-^e". 

La mettiode ICAR [19] utilise quant d elle unlquement les statistiques d'ordre 
4 pour Identifier la matrice G = [gi ... gid des signatures. 
20 En rdsumS, I1d6e mise en ceuvre dans le procddd selon {'invention 

est de construire une observation spatio-temporelle dont les sources 
mSlang^es sont des trains de symboles de I'Smetteur. Chaque train de 
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symboles est par exemple le meme train de symboles decald d'un nombre 
entier de pSriode symbole T. 

Le precede d^crit ci-apres comporte plusieurs variantes de 
realisation dont oertaines sent expliqu6es d titre iliustratif et nuliement 
5 limltatlf. 

Premiere variante de realisation du proc6d§ 

La figure 7 reprSsente un premier exemple de variante de 
rdalisation du pmo666 oO ie signal est regu en bande de base. 

Le precede comporte una 6tape L1 de determination du temps 
10 ^mbole Te en appliquant par exemple un algorithme de detection cyclique, 
tel que celui decrit par exemple dans [1] [1 0]. 

Uetape suivante 1.2 consiste d interpoler les observations x(t) e / 
echantillons par symbole. tel que T=/Te. 

Dans ces conditions oQ fo=0 et bk > Texpression (11) du vecteur devient : 

t(16) 
h(n(l)ITo+jTe)am^(i) + b(mlTe+jTo) pour 0<j<l 

15 Comcne I'dquation (16) est v^rifi^e pour 0 < j < /, le proo^cl^ construit 
I'oiTservation spatio-temporelle suivante (6tape 1.3) d partir des observations 
x(lcTe): 

07) 

hz(n(l)) €Wn<o + bz(ml Te) oD hz(n)= 



2(ml Ta)= 



avec h^j = h(nl To+ jTo) et fe(ml Te)=[ b(ml Te)^ ... b(ml To+( l-1)Te)Y- 
Sachant que x(t) est de dimensiOR-Nx-1 ,-!e vecteur z(t) est alors de dimension 
20 Nb(1. 

Ii(lc) est un vecteur dont la oomposante est le coeffident du fiitre 
filtrant iiniairement le train de symbole {am} sur le n'^"'^ capteur. Le fiitre de 
coefficient vectoriel h(k) depend a la Ibis du fiitre de mise en forme et du 
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canal de propagation. 

Afin d'extraire les U trains de symboles { am-i } d'interet (nombre 
de symboles qui participent d TIES), le precede echantillonne le signal re^u § 
l=(2Lo+1)> en supposant que P < N. 
6 Sachant que le filtre NRZ verlfie 2l-o+1=1 et le filtre de Nyqulst 

2Lo+1=3 pour un roll-off de 0.25, les trains de symboles peuvent etre extraits 
pour oes deux filtres de mise en forme respectivement lorsque P < Nl et 3P 
<NI . 

Le vecteur d'observations z(t) ^tant determine, le prao6dd 
10 applique une m^thode de type ICA pour estimer les U trains de symboles 
{ am-i } associds aux vecteurs de canal ^. =6^ (kj ) . 

La j^^'"® sortie des methodes ICA donne le train de symboles {am. j } 
associe au vecteur de canal h^^. Les trains de symboles {a^. j } estim6s 
arrivent dans un ordre different de celui des trains { am-i } en verifiant : 

a„,j=pexpOai)am^ et Kj^^^i}) ^""^^ 

15 Les trains de symboles { a^. j } sont estjmes avec la memo amplitude car les 
trains de symboles {am.i} sont tous de meme puissance en verifiant : 

E[|am-n(1)fi=-..=E[|am-n(Lc)ft- 

Uetape suivante 1.4 du precede a pour objectif d'ordonner les Lc 
sorties (dm, j. Kj ) dens le meme ordre que les entries (am-i « hz(i)) afin 
20 d'obtenir les vecteurs de canal h^j=h^(kj). Pour oela, le proo6d6 intercorrdle 
deux a deux les sorties am,i et amj en calculant le critgre Cij(k) suivant : 

E[^^,^«>..,*J (19) 



cu{k) = 
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Lorsque la fonction |cij(k)| est maximum en k=kmax les i'**"^ et f sorties 
vdrifient: §m.F SinHanaxj- L'algorithme de classement des sorties am.n(i) ... 
Sm.n(u) est par exemple compose des etapes suivantes : 
Etape n °A.1 : Determination de la sortie dm.imax associSe au vecteur de canal 

5 hj^ioax plus fort module. 

Etape n °A.2 : Pour toutes les sorties &m*j oil j^imax determination des indices 
k=kj maximisant le critSre | cimaxj(k)|. On en dSduit pour cheque j que Sm. imax= 
am-k|.j. Sachant que Cimaxj(ki)= expGaimax- joj) on remet la j**"® sortie d la m§me 
pliase que la Imax'^™ sortie en effectuant am- 14 = cimaxj(l<i) am.j- On remet aussi 

10 en phase les vecteurs de canal en effectuant : h^ikj)= h^ j Cimaxj(i^) • 

Etape n ?A.3 : Cette 6tape remet en ordre les sorties am- 14 et les vecteurs de 
canal h^j^h^Qcj) dans Tordre croissant des Iq sachant que dm = Sm ,imax et Que 

En sortie de ces trols Stapes, on obtient les trains de symboles 
15 { am- k } assocl6s aux vecteurs de canal h,(&p. Sachant que les symboles 
estimds v§rifient am- k= expQ a imax) am-k, la demidre 6tape du procSde va 
consister a estimer cette phase aimax. Pour cela on identifie tout d'abord la 
constellation des symboles ak parmi une base de donn§es composee de 
I'ensemble des constellations possibles. Cette base est constitute des 
20 constellations connues tel que nPSK, n-QAI^. A chaque fois que Ton 
detectera ou que Von aura connaissance d'une nouvelle constellation on 
enrichira la base. 

La figure 6 represente un exemple de constellation de 8-QAM 
lorsque aimax=0 et aimax^. Dans cet exemple de mise en oeuvre, le proctde 
25 comporte alors les Stapes suivantes : 

L'Stape 1.5 suivante consists d ddterminer la phase de sortie 
associde au vecteur de canal de plus fort module. Pour identifier la 
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constellation et determiner la phase le procede effectue par exempie les 

etapes suivantes : 

Etape 1.5 = Etapes B.1, B.2 et B.3 

Etape n'B.I : Estimation des positions des etats de la constellation (points 
5 rouges sur la figure) par la recherche des maximums de Thistogramme 2D 
des points Mk=( reel(dk}> imag(§ic) )• Pour une constellation d M 6tats on 
obtient M couples (Qm, v„) pour 1<m<M. 

Etape n'B.a: Determination du lype de la constellation en comparant la 
position des etats (Qnn v„) de la constellation des {die} a une base de donnees 

10 composee de I'ensemble des constellations possibles. La constellation la 
plus proche est composee des etats (Um.Vm) pour 1<m<l\/l. 
Etape n'B.S : Determination de la phase aimax en minimisant au sens des 
moindres carres le systeme d'equations suivant : 

Qm = COS(<X|max) Um " Sin(afmax) Vm Ot v„, = Sin(ainiax) Um + COS(aimax) Vm 

15 pour 1<m<M 

Le procede peut comporter une etape d'estimation des 
parametres du canal de propagation eh angle Bp et retard xp de Pequation (8) 
par I'algorithme propose dans [8]. Uetape oonsiste k extraire dans un 
premier temps les vecteurs h(n[ Te+ JTe) pour 0<j<l des vecteurs de canal 

20 des vecteurs h^(nj) defini dans {'equation (17). Puis de construire la 

matrice H=[ h(n(1)l Te)..- h(n(Lc) I Te)] de (11) avec les h(nl Te+ jTe) pour 
appliquer la methode [8] d'estimation parametrlque des multi-trajets : (Op, xp) 
1<p<P- 

Deuxieme variante de realisation du precede 
25 Les figures 8 et 9 schematisent une autre variante de realisation 

pouvant comporter deux variantes correspondent respectivement au cas des 
multi-trajets decorreies et au cas des multi-trajets correies par groupe. 
Cas des multi-trajets d6corr6ies. 
Le signal est re^u en bande de base avec { b^i } 
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Les multl-trajets dont les retards verifient | xi^xi | > (2Lo+1) T, ont 
Tavantage d'Stre d§corr6l6s entre eux en vdrifiant : E[s(t-^i) s(t-Tjn=0. En 
observant T^quation (4), on constate alors qu'il suffit d'appliquer une 
mdthode de type ICA iorsque P<N sur I'observation x(t) pour obtenir les 
signaux s(t-a:p) de chacun des multi-trajets. Apres Testimation des signaux 
des diff§rents multi-trajets, le precede determine leurs puissances pour 
garder le signal s(t^pmax) du multl-trajet de plus forte amplitude ppmax- Pour 
determiner ce trajet principal on utilise le fait que les sorties des mdthodes 
ICA verifient en asymptotique : 

sO-Tj ^ (20) 



x(t)=X Ppa(ep)s(t-cp) = X ^i^i(0 avec 5,(r) = 



et 



ai=^Pp a(ep) 

10 ou Yp=Pp^ E[|s(t-^p)|^]. Comme les vecteurs a(ep) sent normes en v6rifiant 
a(Bp)" a(8p)=N, le trajet d'amplitude maximum sera associd a la 

^msB^""^ SOrtiG 

oO aima)r= Siimax" §bnax est maximum. Comme d'apr&s l'6quation (3) la sortie 
*im« (0 =s(t-xpmax) v6rifie : 



^j-axC"™ 'Te+jTo)= ^ hR)(nlTe+jTs-rp™ax)a„« telque 0<j<l 



(21) 



on peut constituer le vecteur observation suivant : 



z(mi Te)= 



L5,^(mn;+(/-i)r,)j 



= il hz(n)amH» ou th(n)= 



(22) 



15 oO hnj=liFo(n I Te+ jTe-tpmax). D'apr&s le module de I'^quation (22), il suffit 
d'appliquer une m6thode de type ICA sur I'observatlon z(ml Te) pour estimer 
les 2Lo+1 trains de symboles { anwi } avec - l-o<n<Lo. Pour extraire les 
incidences Op du canal de propagation, il suffit d'aprds (20) de chercher pour 
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chaque signature ai (1<i<P) le maximum du critere c(G)=| a(8)" §1 Pour 
extraire les retards x\ ^1 du canal de propagation, il suffit d'apres (20) de 
chercher pour chaque signal 5,(0 (1<i<P) le maximum du crit§re 

c(x)=|^,(f-r) 5,(0*1''. 
5 En rdsumd oette varlante comporte par exemple les dtapes 

suivantes : 

Etape n''ll.a.1 : Determination du temps symbole T en appliquant un 
aigorithme de detection cyclique comme dans [1][10]. 

Etape n°ll.a.2 : Echantillonnage des observations x(t) a / echantillons par 
10 symbole tel que T=/ T©. 

Etape n ""ll.a.S : Application d'une m^tiiode ICA sur les observations x(t) pour 
obtenir S,(0 et ai pour 1<i<P- 

Etape n°ILa.4 : Detemiination de la sortie 1= imax ou ai= a" ai maximum. 
Etape n^lLa-S : Constitution du vecteur observation z(t) de (22) d partir du 
15 signal S.^it). 

Etape n'^il.a.e : Application d'une mdthode ICA pour estimer les trains de 
symboles { anvn } oQ -l-o<n< U- On choisit parmi les trains de symboles celui 
qui est associe au vecteur hz(n) de plus fort module : {dm}. 
Etape n^'ILa.?: Determination de la phase Oimax de la sortie associde au 
20 vecteur hz(n) de plus fort module en appliquant les stapes B.I , B.2 et B.3. 
Etape n^'ll.a.S : Remise en phase du train de symboles {§m} en efFectuant 

5„= §m exp(-jaimax) - Le train de symboles constitue la sortie du 

demodulateur de ce sous-proced§. 

Etape n ""li.aS : Estimation des paramdtres du canal de propagation en angle 
25 Op et retard Xp en maximisant pour 1<i<P les critdres | a(e)" ai et 
\Si(t-T) S^iO'*' f pour respectivement les angles et les retards. 
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Gas general avec des multi-trajets quelconques ou corr6les par groupe 
Dans cette variante dont le schema est donne en figure 9, le 
procddd considdre qu'une partie des multi-trajets sont correles. En 
considdrant que I'dmetteur est regu suivant Q groupes de multi-trajets 
corrdlds, le vecteur signal re^u par les capteurs de T^uation (4) devient : 

x(t)=i: i: pp.qa(ep.q)s(t.Tp.c,)+b(t)=i; A,,i^s(t.T ^^^^ 



Oa Aq=[ a(ei.q)... a(epq.q)], Qq=diag([p1^ ...p Pq.J) et S(t, lq)=[ S(t-Tl.q)... s(t- 

"i^Wr avec 2cr[Ti,q ...Tpq,J^. En appliquant une mdthode ICA on estime en 
sortie du s6parateur les signaux et les signatures suivantes d'apr&s : 

(24) 

A=[§i... §pq.q1=[ Ai Ui ... Aq Ucj n et 8(o=n 



oD n est une matrice de penmutation, Uq Vq =aq et Vq E[s(t, t q) s(t, x q)"] Vq" 
=Ipq. Ainsi les trajets d^correles tel que E[s(t-^p^) s(t-Tp',q*n=0 sont regus sur 
des voies ^^(0 et Sj(t) diffdrentes. Les trajets oonrelds ou E[s(t-xp,q) s(t- 

xp\q)*]^0 sont melanges sur une meme vole (t) et sont presents sur Pq ^ la 

fois. Dans la l^^'^ etape de ce sous-procede nous utilisons ce r^ultat pour 
identifier les Q groupes de multi-trajets corrSI^. En prenant les sorties I et J 
du separateur, les deux hypotheses suivantes peuvent etre testees : 



Ho 



[sjiO^bjit) 



Hi: 



(25) 



20 



OU E[bi(t) bj(t-a)*]=0 quelque soit la valeur de t. Ainsi dans Thypothese Ho il 
n'existe pas de multi-trajets comnrums aux deux sorties I et j et dans 
Thypothdse Hi il y en a au moins un. Le test va consister k determiner si les 
sorties s^it) et Sjit-r) sont correl§es pour au moins une des valeurs de x 

v6rifiant | x | < xmax- Pour cela on peut appliquer le test de Gardner [3] qui 
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compare d un seuil le rapport de vraisemblance suivant : 




(26) 



Vij(x)=-2K ln(1-- 



r«(0)r,(0) 



OO Vg(T)<ii => hypothdse Ho 
EtV|j(T)>T| => hypothdse Hi 

Le seuil r\ est d^ermfne dans [3] par rapport d une loi du chi-2 d 2 degr^s de 
liberie. On cherchera tout d'abord ies sorties assoctdes avec la 1^^^® sortie en 
langant le test pour 2<j<PcpcQ et i=1. Puis de la liste des sorties on enldvera 
5 toutes oelles associ^es Si la qui constituera le 1^*" groupe avec q^^L On 
reoommencera la mdme sdrie de tests avec Ies autres sorties non corr6i§es 
avec la 1^ sortie pour oonstituer le 2'^"^® groupe. On effectuera cette 
operation jusqu'au dernier groupe oD au final il ne restera plus aucune vole 
de sortie. On obtiendra finalement en sortie du tri : 



en appliquant ralgorithme MUSIC [1] sur la matrice Aq Aq*^. De ces 
goniometries on deduit Ies matrices Aq. Sachant que Xq(t)= Aq s^(r)= Aq Qq 

s(t, X q), on en ddduit s(t, x q) d une matrloe diagonate prSs en effectuant 
s(t,xq)= Aq%(t). Comme Ies dl^ments des 3(t^) sent composes des signaux 
15 sCt-x p',q), on ddtermine Ies reterds Tp,<r^i.i en maximisant Ies criitdres 
^(-^ )=\Kp(f-^^ Ki(f)* f Kp(f) ^ '21 P*^^ composante de 8(t^q). 
Sachant que E[s, (0 (r) "l=lpq, que Aq "Aq = N Ipq et que Aq (r) = Aq Qq s(t, 

X q), on en ddduit que le groupe de multi-trajets associSs aux amplitudes Qq 
Ies plus importantes maximise le critere suivant : cri(q)==traGe(Aq" Aq). On en 
20 d^ult alors la meilleure sortie assodee a Aqmax et S^^^Cr). Comme d'apr&s 



Aq=AqUq Ot (0 = Vq S(t. X q) pOUf (1 <q<Q) 



(27) 



10 



Les incidences ep,q sont determinSes d partir des Aq pour (1<q<Q) 
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r^uation (3) le vecteurs(t, Tqmax) v^rifie : 



S(m I To+ jTe, Tqmax)= 



pour 



j(mzr.+/r.-r,^ ) 



(28) 



— ^ hpoCn I To+ jTo, 1 qmax) anwi 



0 < j < I et aveC hFo(n 1 Te+ jTe, X qmax)= 



on peut constituer le vecteur observation suivant d'aprds (27): 



z(ml To)= 



s,^(»i/r.+(/-i)rj 



= i) hz(n)aiiwi oQ hz(n)= 



(29) 



ou hnj=Vqmax hFo(n I Te+ jTe,Tqmax). D'apres le modele de TSquation (29), il 

suffit d'appliquer une methode de lype iCA sur Tobservation z(ml Te) pour 
5 estimer les 2Lq+^ trains de symbole { anwi } tel que - LQ<n<Lo. 

En r^umd oette variante comporte les stapes suivantes : 

Etape n''ILb.1 : Ddterminafion du temps symbole T en appllquant un 

algorithme de detection (^clique comme dans [1][10]. 

Etape nH.b.2: Echantillonnage des observations x(t) d I ^chantillonsi par 
10 symbole te! que T=l T©. 

Etape n^'ll.b.S : Application d'une mdtliode ICA sur les observations x(t) pour 

obtenir s(t) et A de l'6quation (24). 

Etape n*'ILb.4 : Tri des sorties suivant Q groupes de multi-trajets correles 
pour obtenir Aq et (0 pour (1 <q<Q) : Pour cela test de correlation de tous 

15 les couples de sorties i et j par le test d deux hypotheses de r^quation (26). 
On cherchera tout d'abord les sorties associ^es avec la 1^^ sortie en langant 
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le test pour 2<j<PQxQ et i=1 . Puis de la liste des sorties on enlevera toutes 
celles associees a la V^'^ qui constituera le 1^®*^ groupe avec q=1. On 
recommencera la meme serie de tests avec les autres sorties non corrdlees 
avec la sortie pour constituer le 2^^"^® groupe. On effectuera cette 
5 operation Jusqu'au dernier groupe ou en final 11 ne restera plus aucune voie 
de sortie. 

Etape n^'ILb.S: Ddtermlnation du meilleur groupe de multl-trajets oQ 
trace(Aq" Aq) est maximum en q= qmax. 

Etape n''ILb.6 : Constitution du vecteur observation z(t) de (29) d partir du 
10 signal s^^CO- 

Etape n'll.b.y : Application d'une mSthode ICA pour estimer les trains de 
symt>oles {am-n} ou -Lo<n< U- On choisit parmi les trains de symboles celui 
qui est associd au vecteur hz(i) de plus fort module : {am-t}. 
Etape n^'ll.b.S : Determination de la phase aimax de la sortie associee au 
15 vecteur hz(i) de plus fort module en appliquant les Stapes B.1 , B.2 et B.3. 
Etape n^'ILb.S : Remise en phase du train de symboles {dm} en effectuant 

dm axp(-jaimax) - Le train de symboles {&^} constltue la sortie du 
demodulateurde ce sous-procede. 

Etape nll.b.10: Estimation des parametres du canal de propagation en 
20 angle Oq,p et reterd Tq,p. Les incidences Op^q sent ddferminSes d partir des Aq 
pour (1<q<Q) en appliquant I'algorithme MUSIC[1] sur la matrice Aq Aq" De 
ces goniometries on en dSduit les matrices Aq pour en ddduire une estimd de 
s(t> Tq) en effectuant S(t.Tq)= Aq^Xq(t). Comma les 6ldments des s(t,Tq) sent 
composes des signaux s(t-Tp.q), on determine les retards xp.q-xi.i en 
26 maximisant les critdres c(t )=| j^CO*!^ ou ^^.p(0 est la p'^™ 

composante de s (t,Tq). 

Autre variante de mise en osuvre du proc6d6 

Estimation de la frequence porteuse et deduction des {am}- 

Cette technique va consister d estimer la frequence porteuse fo de 
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i'^metteur ou le complexe zo= exp(j27tfoTe) pour ensuite deduire les symboles 
{am} des symboles {bm} en effectuant d'apr§s (3) : 

am=bm exp(i27tfomlTe)= bm zo"^' (30) 

Cette etape est appliquee apr^s T^tape n''l.4 de remise en ordre des 
symboles et des vecteurs de canal. D'apr^s les equations (3) (17) et (7) (8), 
5 on dispose des vecteurs de canal suivant : 

(31) 

h^(n) = *^ -v'.»*e-*ey pourneQ 



«/+(/-!) 



h(nn;+(/-i)rj 



oD CI ={lnd[q>4.o. rp4Lqi(n)=1 pour un p tel que 1<p<P} 

Sachant que CI ={ n^ <...<nix^, on constitue d partir des vecteurs h,(n)un 
grand vecteur b tel que : 

(32) 



w = 



h,(«i) 

KM 



.KM. 

La recherche de fo va consister a maximiser le crlt^e suivant : 
Porteuse(fo)=|w" c(exp027cfoTe)}|^ 



(33) 



Avec c(zb)= 



et oCi c(n,zo) 



ni 



2o 



nl+l 



_Zo 



«/+(/-!) 



10 Les stapes du proc6d§ adapts au cas d'un dmetteur k lT§quenoe non nuile 
sont les sulvantes : 
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Etape n°lll.a.1 : Etapes 1.1 jusqu'ii i.4 decrites ci-dessus pour obtenir les 

trains de symboles {b„_t^ } associds aux vecteurs de canal h,(ftp . 

Et&pe n°lll.a2: Constmction du vecteur w de I'iquation (32) § partir des 

5 Etape nlil.a.3 : Maximisation du crit^re Porteuse(fo) de l'6quation (33) pour 
obtenir fo. 

Etape n1ll.a.4 : Application de {'equation (30) pour d^duire les symboles {am} 
des symboles {bm}. 

Etape n*1ll.a.5 : Etapes I.Sjusqu'a 1.7 precedemment decrites. 
10 Dans ie cas d'un ^metteur d frequence non nulle et pour des multi- 

trajets decorrSlte les Stapes sont les sulvantes : 

Etape n*1ll.b.1 : Etapes ll.a.1 Jusqu'd li.a.4 decrites ci-dessus pour obtenir Ie 
vecteur z(t) de I'Squation (22). 

Etape n°lll.b.2 : Application des mSthodes ICA [4] [5] [15] [19] pour estimer 
IS les Lc trains de symboles {b„ j} associ6s aux vecteurs de canal . 

Etape n^lll.b.S: Remise en ordre des trains de symboles {b^j} et des 
vecteurs de canal en appliquant les Stapes A.1. A.2 et A.3 pour obtenir 
les trains de symboles {b^t^ } associds aux vecteurs de canal =b,(*y) - 
Etape n°ili.b.4 : Construction du vecteur w de I'Squation (32) ^ partir des 

20 h,(fc^). 

Etape n'llLb.S : Maximisation du critSre Porteuse(fo) de I'Squation (33) pour 
obtenir fo. 

Etape n°ill.b.6: Application de i'equation (30) pour deduire les symboles 
{am} des symboles {bm}- 
25 Etape n <>lll.b.7 : Choix du train de symboles assodS au vecteur hz(i) de plus 
fort module : {S^mi- 

Etape nlll.b.8 : Etapes li.a.7 jusqu'd ll.a.9 dScrites pr6cSdemment. 

Dans Ie cas d'un Smetteur d frequence non nulle et pour des multi-trajets 
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corrdlds les etapes sont par exemple les suivantes : 

Etape n**IILc.1 : Etapes ILb.1 jusqu'a ll.b.6 n*'Z2 pour obtenir le vecteur z(t) 
de rdquation (29). 

Etape : Application des m6thodes ICA [4] [51 [151 [19] pour estlmer 

s les U trains de symboles {b^j} assocles aux vecteurs de canal h,^ . 

Etape n°lll.c.3: Remise en ordre des trains de symboles {b^j} et des 
vecteurs de canal h^ j en appllquant les stapes A.1 , A.2 et A.3 afin d'obtenir 
les trains de symboles {4,-*^ } associe aux vecteurs de canal h,j=hz(*/)- 
Etape n°IILc.4: Construction du vecteur w de Tequation (32) d partir des 
10 h^Qcj). 

Etape n*1ILc.5 : Maximisation du critdre Porteuse(fo) de I'equatlon (33) pour 
obtenir fo. 

Etape n'^lll.c.e : Application de I'dquation (30) pourd^duire les symboles {am> 
des symboles {bm}. 

15 Etape n1ll.c.7 : Choix parmi les trains de symboles celui qui est associ^ au 
vecteur hz(i) de plus fort module : {am-J- 

Etape nlll.c.8 : Etapes ll.b.8 jusqu'd ll.b.10 pr§c6demment ddcrites. 
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REVENDICATIONS 



1 - Proc6d6 de demodulation aveugle d'une source ou 6metteur de forme 
5 d'onde Hn6aire dans un systems connportant une ou plusieurs sources et un 

r^seau de capteurs et un canal de propagation caract^ris^ en ce qu'li 
comporte au moins les stapes suivantes : 

• d6terminer le temps symbole T et on echantiilonne d Te tel que T=ITe 
(I entier), 

10 • d partir des observations x(kTe}, construire une observation spatio- 
temporelle z(t) dont les sources m^iangSes sent des trains de 
symbole de Tdmetteur , 

• appliquer une mdthode de type ICA sur le vecteur d'observation z(t) 
pour estlmer les U trains de symboles { SmA } associ^s aux vecteurs 

15 de canal h. y = h^ikj) , 

• ordonner les U sorties (am. j, fi^j ) dans le mSme ordre que les entries 
(am-i , hz(i)) afin d'obtenir les vecteurs de canal de propagation 

• dSterminer la phase aimax associSe aux sorties. 

20 

2 - Precede selon la revendication 1 caract6ris6 en ce que Ton estime les 
parametres du canal de propagation pour determiner la fr§quence porteuse 
afin de compenser les trains de symboles pour les obtenir en bande de base. 



25 3 - Pmc6d6 selon la revendication 1 caractdrisd en ce qu'il oomporte une 
§tape d'estimation des param&tres du canal de propagation en angle Op et 
retard tp,. 
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